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1 Introduction 



The Gross-Neveu model |IJ, is an example of a fermionic field theory which ex- 
hibits a non trivial ultraviolet (UV) fixed point in 3 dimensions. From a field 
theoretic point of view such models, the non linear sigma-model is another one, 
provide at the fixed point, examples of quantum field theories (QFT) which are 
interacting and scale invariant. In the Euclidean formalism such models describe 
classical statistical systems undergoing second or higher order phase transitions. 
Euclidean scale invariant QFT are then relevant from a phenomenological point of 
view, since they describe experimentally accessible phase transitions phenomena 
and the corresponding critical exponents of such transitions have been calculated 
with great accuracy in many realistic cases. With this respect, results have been 
obtained mostly for theories which have a larger symmetry, conformal invariance 
(particularly in 2 dimensions where the conformal algebra is infinite dimensional). 
This is also the case for the nonlinear a model in 3 dimensions, in the large N limit, 
as was shown in Refs.|2|, BL We will see instead, that the 3 dimensional GN model 
in the large N limit is not conformally invariant at the fixed point. 

In this paper we study the GN model, on 3-d manifolds of constant, non-zero 
curvature. The model exhibits on R 3 a two-phase structure, the phase transition 
occurring for a non-zero value of the coupling constant (see for example Q for 
a review). At the critical point there is a symmetry breaking and the fermions 
become massive. The model furnishes a description of phase transitions in classical 
superconductors, as it was shown first in ||, where the stability conditions for the 
effective action in the large N limit are seen to imply the BCS gap equation. 

We consider the model on manifolds of the form M. = S x Si and analyze 
the limit Sl—*R as /3 goes to infinity; £ is a 2— d manifold with constant, non- 
zero, curvature. The interest for Euclidean field theories living on manifolds of 
the form M. = £ x Si is due to the fact that the radius of the circle Sp can be 
interpreted as the inverse temperature of some two-dimensional statistical system. 
The action describing the model is then dependent on the coupling constant present 
in the theory and on the parameter (3. Phase transitions can occur with respect 
to both the parameters. There is also a phenomenological reason to study critical 
phenomena in curved spaces. In fact, applying to a system an external stress it is 
possible to deform its microscopic structure in such a way to change the effective 
distance between points. At a critical point, universality suggests that the details 
of the microscopic structures do not matter; but the system is still sensitive to the 
deformation through the effective metric tensor density. Nevertheless, the study 
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of quantum effects for field theory at finite temperature in curved space is an 
interesting subject by itself, since only in few cases the quantum theory can be 
computed and in perspective one could imagine some applications to cosmological 
scenarios. 

There exists a very reach literature on the 3-d GN model. The issue of critical 



exponents and (3 function is addressed for example in Refs.0-|jTO|- In Refs. UTTj] |[T3|| 
the effects of an external electromagnetic field are considered. The thermodynami- 
cal behaviour and the proof ofl/N renormalizability can be found in Refs. flTl|-[|16| . 



With respect to curvature induced symmetry restoration, see for example Ref.fPT| 
(GN model in 2 dimensions), and Ref.flTB" 



The paper is organized as follows. In section 2 the properties of the GN model 
are reviewed, whereas in section 3 the large N limit approximation is discussed. 
Section 4 is devoted to the flat space analysis and in section 5 we introduce the £- 
function regularization. Thus in sections 6 and 7 the spaces S% x Si and H 2 r x Si are 
analyzed at the critical coupling and the limit (3 — > is discussed. The asymptotic 
behaviour of the correlation function for the above manifolds is studied in section 
8. Finally in section 9 we give some concluding remarks. Some technical details 
are left for the appendices. 



2 The Gross— Neveu model in three dimensions 

In this section we study the Gross-Neveu model on a Riemannian manifold (A4, g). 
The model is described in terms of a 0(n) symmetric action for a set of N massless 
Dirac fermions. The Euclidean partition function of the GN model in 3-dimensions 
in the presence of a background metric g^ v {x) is given by 

Z[g] = J V[$\ V[$\ exp {- d 3 x ^g ^(xW^x) + |(^) 2 ] } , (2.1) 

where % — 1, 2, • • • , N, yis the Dirac operator on A4, and q is the coupling constant^. 

The Dirac matrices on a generic 3-d manifold Ai are given in terms of the Pauli 
matrices a a by the expression 

7 M = V^ a a a , with fj,, a = 1, 2, 3 (2.2) 

where V^ a denote the dreibein defined by the equation 

5> = V^ a (x)V Utb (x)5 ab . (2.3) 

*According to our notation the Dirac matrices obey the following algebra: 7 M = 7L {7^, 7,,} = 
2g^, and Tr (7^) = 0. Thus, the Dirac operator is antihermitian = — y. 
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The covariant derivative V„ acting on a spinor field is denned as ML [19 



v„ = <9 M + r M (x) 



where is the spin connection 



(2.4) 



(2.5) 



The kinetic term of the action (|2.1|) is invariant under the conformal transformation 
of the metric, g^ u {.x) — > fl 2 (x)g fiu (x). This can be easily checked observing that the 
Dirac operator transforms according to [IS] 



f[Q(x)] 



and assuming the spinors to be conformal densities of weight —1 

ipi(x) -> Q' 1 (x)ipi(x) . 



(2.6) 



(2.7) 



On the contrary, the interacting term in (|2.1| ) violates conformal invariance. 



An equivalent way of rewriting the partition function ( |2.1|) , but more suitable 
for our purposes is obtained by introducing an auxiliary scalar field a, such that 



Z[g] = J V[$\ Z#] V[a] exp |- d 3 x 



2q 



-a 



(2.8) 



The new field has no effect on the dynamics of the theory, since from the point of 
view of functional integration the integral over a merely multiplies the generating 
functional by an irrelevant constant. We regularize the generating functional Z in 
the ultraviolet by introducing a cut-off, A, in the momentum space. Before doing 
that, let us compute the canonical dimensions of the fields and of the coupling in 
our action. In mass units they result to be 



1 



01 



1 



(2.9) 



By replacing the dimensional coupling constant l/g(A) with the dimensionless ratio 
A/g(A), the regularized partition function can be formally rewritten as 



Z[g,A}= / V A [ij] VM V A [a] exp 



M 



d 3 x ^fg 



(2.10) 



where Pa^] — II dip(k) and similarly for the other fields. 

|fe|<A 



4 



3 The large N limit 



The GN model is not exactly solvable, so that we have to use some approxima- 
tion method. As for the non linear cr-model, the existence of a non trivial UV 
fixed point shows that the large momentum behaviour is not given by perturba- 
tion theory above 2 dimensions, where the theory is asymptotically free (see for 
example Ref.Q). Other techniques are required, like the 2 + e expansion, which 
relies on the fact that the 2-d model is renormalizable in perturbation theory, or 
the 1/N expansion, which we will use in the paper. The model has been proven to 
be renormalizable in 3-d in the 1/N expansion 

In this limit, which means N oo keeping Nq(A) fixed, the generating func- 
tional can be calculated using the saddle point approximation. For this purpose 
we integrate over N — 1 fermion fields, rescale the remaining fields ipN, 4>n to 
\/N — 1 ip N , \/N — 1 ijj N , respectively, and redefine (N — l)q(A) as 5(A). Thus 
we get 

Z[g,A,q(A)} = I V A [i> N ] VSn] V a [<t] exp {-(N - l)Tr log A (y+ a)} 

A 

IM 



x exp \ —(N - 1) / d 3 x ,/g 

JM 



^j N (f+a)ij N - —a 2 (x) 
zq 



(3.1] 



In the limit N 00 the dominating contribution to the functional integral comes 
from the extremals of the action. For an arbitrary metric g^x), these are obtained 
by extremizing the action with respect to ^at(x) keeping a(x) and ipN fixed and 
vice-versa. Hence, a set of equations (gap equations) is obtained 

tp N (f-o-) = , (3.2) 
(f+a)^ N = , (3.3) 

4>n*Pn = -ttto- - G A (x,x;cr,g) , (3.4) 
q(A) 

where G A (x, x\ cr, g) = (x\( J f+ cr) _1 |a;)A is the two-points correlation function of 
the f/>Ar-field, evaluated for x — > x' . Note that the derivative in Eq. (|3.2j ) is acting 
on the left. 

In the following analysis the generating functional Z will be evaluated in the 
large A^ limit, at the uniform saddle point 

(cr) = m , (ip N ) = b , (^jy) = b . (3.5) 

The values m, b, b will be given by constant solutions of gap equations. The 
quantities b and b represent the vacuum expectation value (v.e.v.) of fermion fields, 
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while m, if positive, can be regarded as the mass of the field fluctuations around 
the vacuum. In the language of condensed matter physics b is also addressed as 
spontaneous magnetization. 

Once we find the saddle point solutions of the action, we can compute the 
generating functional of connected Green functions W = — log(Z) (the free energy), 
at the saddle point, to the leading order in the 1/N expansion 



The gap equations are the equations of motion of the classical field theory: the large 
N limit of the Gross-Neveu model. The ground state will then correspond to the 
solution which minimizes the free energy. If the background metric is homogenous, 
we expect the ground state solution of the gap equations to be constant. 



Before solving the gap equations for curved spaces, it is worth noticing that short 
distance divergences of the Green's function G\(x, x; m 2 , g) are independent of the 
curvature of the space (the argument is completely analogous to the one given in 
Ref.fH] for scalar fields). This observation allows us to calculate the critical value 
of the coupling constant in the simple case of flat space R 3 , the critical coupling 
constant being, in fact, the value of q(A) which makes the divergences cancel in 
equation ( |3.4|) . 

In the following we recall the two-phase structure of the GN model in flat space 
using Pauli-Villars regularization. Then we solve the gap equations for the physical 
parameters m, b, and b and use these values (which must be independent of the 
regularization scheme) to determine the critical coupling constant in the (^-function 
regularization. 

The uniform saddle point is determined by the solution of the gap equations 



A f 

W[g, A, q(A)} = N Trlog A (y+m) / Sx^fgm 2 



q Jm 



(3.6) 



4 Flat space (R 3 



) 



mb 








bb 




A 




with 



G\(x, x;m, g) = lim(x|(^ + m) \x') 



Am m 2 



(4.3) 



2tt 2 4tt 
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Substituting the last expression into Eq. (|4.2|) and recalling the definition of critical 
coupling constant we find 



q c 27r 

Posing M = (A/q) - (A/q c ), Eq. $~% ) becomes 



- = ■ (4-4) 



bb = Mm + —m - G A . (4.5) 

Qc 



Equation ( [4.1| ) requires that either m or b be zero, or both. If m is zero, ( |4.2| ) 
implies b = b = 0. If b is zero, we have 

2 

777, 

Mm = . (4.6) 

47T 

Hence, we can conclude that: m and b are zero at the critical point M = 0; 
m/0, 6 = for M < 0, while m = b = when M > 0. The point M = is where 
the phase transition takes place, while for q > q c the fermions acquire mass and 
the symmetry is spontaneously broken. In the next section we use m c = b c = to 
determine the critical coupling constant in the (^-function regularization. 



5 (—function regularization 



We will use throughout the paper the ^-function regularization. In fact, it results 
to be more tractable when we are on curved spaces, even though other regular- 
izations such as the Pauli-Villars have a more immediate physical meaning. Since 
the critical value of the coupling constant at which the theory becomes finite is 
independent of the background metric, we compute this critical coupling on R 3 . 

Let us consider the squared Dirac operator in d dimensions 

K 

T 



A 1/2 ^y=-V? /2 + - . (5.1) 



Note that Ax/2 is not the conformal spin 1/2 Laplacian, which is instead the com- 
bination [22 

where TZ denotes the Ricci scalar. 

Given the eigenvalues, A^ + m 2 , of the operator (Ai/2 + m 2 ) and an orthonormal 
basis of corresponding eigenvectors {ip n (x)}, the local (^function is defined as 

C(s,x)=^l + m 2 )- s \Mx)\ 2 , (5-3) 
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where the sum includes degeneracy, and in case m — 0, X n must be non vanish- 
ing. If the eigenvalues are continuous, the sum is replaced by an integral. The 
Green's function for the operator y + m can be defined through the (^-function 
( |5.3| ). Observing that 

(x\(f+m)- l \x) = Y,( x \^n)M(¥+m)~ l \^i)(Hx) 

n,l 

n X 2 n + m 2 

and recalling that the eigenvalues of the Dirac operator y always appear in pairs 
±iA n , we regularize the Green's function of the operator y + m as 

G s (x, x\ m, g) = m(x|(A 1 / 2 + m 2 )~ s \x) = m((s, x) , (5.5) 

and 

G(x,x;m,g)=m \im((s,x) . (5.6) 

On homogeneous spaces such as the ones we will be considering in this paper, 
((s,x) turns out to be independent of x. 

On R 3 the gap equation ( |3.4| ) becomes, in this regularization scheme, 
1 - f d?k 1 

m ^W) = bb + m ^A s ^ = bb + m ^J (27r )3( fc2+m2 )s 

= 66+^lim r dt^-Jk 2 dke^ k2+m ^ , (5.7) 

Z7T Z s ^ JO 1 [SJ J 

where the regularized coupling A/g(A) in the Pauli-Villars regularization has been 
replaced by l/q(s) in the (^-function regularization. Here we have used the Mellin 
transform to analytically continue the (^-function. Note that 

C(S)= / (2n)%k 2 + m 2 Y ' (5 - 8) 
has no pole at s = 1. It is now easy to verify that 

hm — - = bb - hm » , 2 =bb-\ . (5.9) 

s^iq(s) ( 47r )l r(s) 4vr 

Using the critical values of b, b and ra, found in the previous section (being physical 
quantities they are regularization scheme independent), we get, in the (^-function 
regularization 

— = . (5.10) 
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Hereafter we use this value of the critical coupling, since it is independent of the 
background metric. 

At the critical point the large N limit of the free energy ( |3.6| ) becomes then 

N 

W c ((?) = ylogdet(y+m c ) . (5.11) 

Observing that logdet(y+ m c ) = (1/2) logdet(Ai/ 2 + m 2 ) one can define the free 
energy density w c (g) through the (^-function (|5.3|). We have 



A + / A , 2x I 0. if dim(ker[A 1/2 + m 2 ]) ^ 

det(A 1/2 + m c )H e _ c , (0); ifdim(ker[Ai/2 + mc2]) = • (5-12) 



The free energy density at the critical point is then regularized as 

JO, ifdim(ker[A 1/2 + m 2 ])^0 

Wc[9) ~ I -|C(0), if dim(ker[A 1/2 + m 2 ]) = " 



This quantity is indeed a density, since the ^-function defined in Eq. (|5.3|) contains 
an inverse volume factor through the squared modulus of eigenvectors. 

We wish to stress at this point a major difference with the non-linear sigma 
model. It was found in Ref.0 that the free energy for the non linear a model at 
the critical point is equal to 

W c = — \ogdet(p+m 2 c ) , (5.14) 

where □ is the conformal scalar Laplacian []T9] . Using this fact (of the Laplacian 
being conformal) it was shown that the free energy is conformally invariant. At a 
first sight equation (|5.14|) looks identical to the expression ( |5.13|) . But, we have 
seen that Ai/ 2 is not the conformal Laplacian, then the arguments used in Ref. @ 
do not apply and we can conclude that the free energy of the Gross-Neveu model 
at the critical point is scale invariant but not conformally invariant. 



6 The Manifold S r 2 x S} 

In this section we study the large N limit of the Gross-Neveu model on the manifold 
S 2 x Sa (r and (3 are the two radii). We parametrize this space by x M = (r, x, 9), 
where < r < 2tt, —tt/2 < \ < tt/2, and < 6 < 2n. The metric tensor is then 
defined as 

g^da? ® dx u = r 2 cos 2 x dr ® dr + r 2 d\ ®dx + [3 2 d6 ® dO . (6.1) 
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Using the definition ( |2.2| ), the Dirac matrices on S 2 x Si are given in terms of the 
Pauli matrices a a , by 

7l = rcosx(Ji , l2 = ra 2 , 73 = ^3 • (6.2) 



The spin connection (|2.5| ) results to be 



i 



r^(ar) = --cr 3 sinx 5^ , (6.3) 

while the covariant derivative acting on a spinor field is given by ( |2.4| ). Substi- 
tuting ( |6.2| ) and ( |6.3| ) in the Dirac equation 

y^A = ^A . (6.4) 

we find the eigenvalues of y to be i\f n with 



, (2n + l) 2 (/ + 1) 2 

Kn = ±\ r 4/32 } + ^ = 0,1,- n = 0,±l,.. (6.5) 



and degeneracy 2(/ + 1); note that there are no zero modes. The details of calcu- 
lation are given in Appendix A. 

The eigenvalues of A1/2 + m 2 are then given by |A/ n | 2 + m 2 , with degeneracy 
4(1 + 1). The ^-function (|5.3|) for this operator is 



(M = ^ W z E E 



n=— 00 Z=0 . 



(2n + l) 2 (/ + l) 2 2 
+ ^ + m 2 



(/ + !) , (6.6) 



4(3 2 r 2 
while the gap equations (|3.3|) and ( |3.4j) become 

( 7 % + m)6 = , (6.7) 



66 = mlim<; 



The first one yields 6 = 0, so that we have to solve 



= m lim 



{^)~ C(S,m) } ' (6 ' 9) 



Since lim s ^! l/q(s) = and linx,^! £(s, m = 0) is well defined, m = is a solution. 
For m/0 there is no value of m such that ((s, m) = 0. Then m c = 0, 6 C = 6 C = 
are the critical values of the physical mass and the vacuum expectation value of the 
if}, if) fields on S 2 x Sp, in the large N limit. Note that m c = is not in contradiction 
with the fact that the correlation length has to be finite due to the finite size of 
the background. In fact we have already observed that the operator A1/2 has no 
zero modes, so that the smallest eigenvalue of A1/2 + m 2 is non zero. 



10 



To evaluate the large N limit of the free energy density at criticality, we have 
to take the derivative of the (^-function at zero and substitute the value of m c that 
we have just found, in Eq. ( p,13| ): 



1 f d 00 00 

-20^ x E E 



8n 2 r 2 f3 [ ds 



n=— oo i=0 



(2n + l) 2 (1 + 1 



40 2 



+ 



„2 +m c 



(/ + !). (6.10) 



We first take the Mellin transform of ((s,m 



C(s,m c 



An 2 r 2 T(s) 



r oo 0000 r r / i\ 2 

/ ^ E E^xp - + Z 2 ^ + 



m 2 /3 2 



(6.11) 

then use a specialization of the Poisson sum formula for the sum over I (the formula 
is obtained in Appendix B), which allows us to exchange the sum over I with the 
integral in t 

,3 



E^ ex P H 



i=l 

and we get 



2^0' 



;t~ 3/2 J°° dx x cot(x) exp|-x 2 ^| , (6.12) 



C{s,m c 



8vr 5 / 2 r(s) 



0~ 2s ~ 4 / dtt s ~ 5/2 P dx [xcot(x) -1} 

JO „ J —oo 



x exp 



-2s— 3 />oo 



8vr 2 r(s) 



/ rftt s - 2 ex p 



n + 



+ m 2 /? 2 



/3 2 t_ 

tl=A + B. (6.13) 



Here we extracted the part of the integral in x which would become divergent once 
we exchange the order of integration. We find also convenient not to put m c = 
until the end of the calculation. 

To evaluate the contribution of A to the derivative of the ( function we first 
perform the integral in t. We get 



lim 

s-^0 



ds 



'TO JO X 6 



x {1 + 2x 







n + 



4n 2 r 2 (3 n=0 
1/2 



+ m 2 /3 2 



exp < — 2x 







n + 



+ m 2 /? 2 



1/2' 



(6.14) 

The principal value of the integral can be evaluated using the method of residua, 
since it has only simple poles. Thus, for m c = 0, we find 



lim 



ds 



3Ca(3) 



E 



- m c =0 



32vr 4 /? 3 8vr 2 r 2 /?^ n 



^ -cosech 2 ( nir— 



(6.15) 
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where (r(z) denotes the Riemann (^-function. 

The second integral in ( |6.13| ) is easily computed by using the Poisson formula 
for the sum over n |23| 



J2 ex p < 

n=— oo 

Thus one finds 
B = 

+ 



n 



l + 2£(-l)"exp 



ix 2 n 2 



n=l 



(6.16) 



1 



47r 3 / 2 r(s) 
1 



0- 



(3- 



-2s-3 



00 r°° , ( \ir 2 n 2 

E(-ir/ dtt^ 2 e W \- ^f- 

n=l J ° { I t 



■m 2 c (3 2 t 



m 



3-2s n 



87r 3 / 2 r( s ; 

In the limit s — > m c = 0, we find then 



lim 

s->0 



d_ 

ds 



B 



- m c =0 



— Y 



-D n 



rv 



3Ci?(3) 
32tt 4 /? 3 



(6.17) 



(6.18) 



Summing fl6.15|) with (|6.18|) we finally get the regularized free energy density at 
the critical point 



U'c 



167r 2 r 3 \(3 



^ ) ^— cosech 2 (nir— 



(6.19) 



It can be checked that the series is convergent to a negative value for any finite 
value of (3/r, so that the free energy density is positive definite. The series can be 
evaluated numerically. In Figure 1 we plot the free energy density of Eq. (|6.19|) as 
a function of the ratio j3/r (up to the overall factor l/167r 2 r 3 ) 

The limit j3 — > oo , which corresponds to zero temperature, yields w c = 0, in 
agreement with the result that we would get by direct calculation on S 2 x R. The 
limit r — > oo, which corresponds to the manifold R 2 x S^, yields 

(-1)" = 3^(3) 

n 3 

n=l 



lim w c 

r—>oc 



16tt 4 /? 3 



64vr 4 /3 3 



(6.20) 



where Ch(3) = 1.20205.... This represents the free energy density for the GN 
model on flat space, at finite temperature 1/(3. As we can see it goes to zero with 
temperature. To our knowledge this result was first obtained in [fTI], where the 
free energy density of the GN model is evaluated on R 2 x S* 1 in the presence of a 
magnetic field. The limit of zero magnetic field is in agreement with ( |6.20| ) once 
we observe that m c is zero in this limit. 



Finally, we wish to note that ((0, m c 
of the model at criticality. 



in agreement with the scale invariance 
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7 The Manifold H 2 r x S} 

We consider the product manifold H% x Si. Hf, is a 2-dimensional pseudosphere, 
namely, it is obtained globally embedding a hyperboloid in R 3 endowed with 
Minkowskian metric. Either sheet of the hyperboloid models an infinite space- 
like surface (hence with Riemannian metric) without boundary. This surface has 
constant negative curvature and it is the only simply connected manifold with this 



property P^fl . We parametrize H% as 

E 2 r = {z = (x,y), x e R, < y < 00} , (7.1) 

while the circle S 1 of radius (3 is parametrized as before by 9, < 9 < 2tt. The 
scalar curvature of H 2 is 1Z = — 2/r 2 , where r is a constant positive parameter. 
The metric tensor on the whole manifold is then given by 

gpvdx? ® dx u = T — (dx ® dx + dy g> dy^j + (3 2 d9 ® d9 , (7.2) 

where x^ = (x, y, 9) with /i = 1, 2, 3. The Dirac matrices on H 2 x Sp are given in 
terms of flat Dirac matrices by 

T r 

71 = -ax , 72 = -cr 2 , 73 = (3a-i , (7.3) 
y V 

while the spin connection defined in fl2.5p, is 



r M = ^^ . (7.4) 

To find the spectrum of the Laplacian, we could in principle proceed as in the 
previous case, by using the algebraic method described in Appendix A. However, 
in order to write the ^-function all that we need is the heat kernel of the Laplacian. 
Thus, we use the algebraic method only to establish that the eigenvectors of the 
Dirac operator (and hence of its square) are of the form 

\i>n,k) = \Xk) ® |0n) , (7.5) 

where 

\<f> n ) = exp f z (n + I) 9} (7.6) 



^3 r lV 2 

are scalars depending on 9 only, while \xk) are spinors depending on x, y. They 
are more explicitly described in Appendix A. This means that the Laplacian can 
be factorized as 

A 1/2 (tf f 2 x 5J) = A 1/2 (H?) + Ao(5j) (7.7) 
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and the heat kernel is just the product of the heat kernels of the Laplacians on H 2 
and Sp, respectively, 



h(t;x,y,0) = h H ?(t;x,y)h s i(t;9) 
The heat kernel for the spin 1/2 Laplacian on H 2 is p5| 



h H ?(t; z, z') 



00 w cosh(it;/2) 

dw 



y / cosh(c//r) + 1 ^ d / r y / cosh(w) — cosh(d/r 
where d is the geodesic distance among points on H 2 , 



exp 



— w 2 r 2 
4t 



(7- 



(7.9) 



cosh 1 - 



1 + 



\z — z 



2yy' 



t\2 



and z = x + iy. In the limit d — > Eq.([7.9|) becomes 



hm(t;z = z') 



2(7Tt)3/2 Jo 



du> u> coth(w) exp 



-w 2 r 2 



(7.10) 



(7.11) 



The equal points heat kernel of the scalar Laplacian on Si (the spectrum is found 
in Appendix A) is 



1 00 

hsl ^ e = e ' ) = 2^3 £ ex P 

M n=— co 



n + 



^ 2 



(7.12) 



We are now ready to write the spectral (^-function for the operator Ax/2 + m 2 on 
the product manifold. It is defined in terms of the heat kernel as 



I f°° r -\ 

C{s,m) = f^J o dt t s ^h H 2{t- z = z')h s i{t; 6 = 6') exp {-m 2 t} 

On H 2 x Sp the gap equations ( |3.3|) and ( |3.4|) become 



( 7 T M + m)b 
bb 





mlim< 



The first one yields 6 = 0, so that we have to solve 

= m lim(— — — ((s,m)\ 
s-+i(.q(s) J 

By rescaling t — > t/3 2 , the ^-function ( [7.13|) becomes 



(7.13) 

(7.14) 
(7.15) 

(7.16) 



C(s,m) 



4vr 5/2 r (< 

00 

£ exp 



/■CO [ /"C 



dw w coth(io) 



w 2 r 2 



- I (n+~) +m 2 l3 2 \t 



(7.17) 
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One can check that, since lim s ^ C( s ; m — 0) is finite, m = is a solution of (|7.16| ). 
Thus, at the critical point we have 







(7.18) 



To find the free energy density at criticality we have to consider the derivative 
of the (^-function at s — 0, and then evaluate it for the critical value of m just 
found. The calculations go along the same lines as those for S 2 X Sh each time we 
exchange integrals and series we have to verify that no divergences are introduced, 
and in that case eventually regularize them. 



We have 
C(s,m c ) 



4vr 5 /2r(s 
x exp 



2s-4 



dt t s ~^ 2 ^2 dw [w coth(w) - 1] 



w 2 r 2 



(3 2 t 



n+~) + m 2 c /3 2 



POO °° POO 

+ / dt t s - 5/2 Y / dw exp 
= A + B . 



w 2 r 2 



(3 2 t 



n + 



+ m 2 c (3 2 



(7.19) 



The derivative of the first integral in ( [7.19| ), evaluated at s — > 0,m c = 0, gives 



lim 

s->0 



d 
ds 



- m c =0 



^7r 2 r 3 Jo w 



°° dw 

2 



coth(w) 

w 



wr , ( wr\ 
— coth — 







y 



sinh (w r/(3) 
(7.20) 

The derivative of the second integral in (|7.19|) , at s — > 0, m c = 0, gives 



lim 

s^0 



ds 



m c =0 



E 

n=l 



(7.21) 



Putting together ( |7.20 ) and ( 7.21 ) we get the expression of the free energy density 
at the critical point to be 



1 



W f . 



16n 2 r 3 [ 
3Cfi(3) (r 



00 dw 
w 2 

3 



coth(w) 

w. 



wr 



wr 



— — coth — — 







r//3 



47T S 



sinh (w r/0) 

(7.22) 



The integral can be evaluated numerically and the free energy density w c results 
to be positive definite. In Figure 2 the free energy density of ( |7.22| ) is plotted as a 
function of the ratio P/r, up to the overall factor l/167r 2 r 3 ). 
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As for the manifold S 2 x Sp, the limit of zero temperature can be evaluated 
analytically and we get 

lim w c = . (7.23) 

The limit r — > oo, which corresponds to the manifold R 2 x S 1 yields 

r 3Cr(3) 

hm w c = - . „ . 7.24 
As expected, we get the same result that we found in (|6.20|) for the limit S 2 — > R 2 . 



8 Asymptotic behaviour of the correlation func- 
tion 

To understand the effects of the curvature on the second order phase transition 
exhibited from the model on flat space, we investigate the behaviour of the two- 
points critical correlation function as the distance among points goes to infinity. 
The correlation length which characterizes such a behaviour, can also be regarded 
as the inverse of the smallest eigenvalue of the operator under consideration. For 
this reason it is then equivalent, but easier, to study the correlation function of the 
squared Dirac operator instead than the Dirac operator itself. 

When a second order phase transition occurs, the correlation length diverges 
and the two point Green's function follows a power law for large distances. This 
is what happens for the GN model on R 3 at the critical point. We don't analyze 
the case of S 2 x Si, because finite size of the manifold in all directions forces the 
correlation length to be finite. On the other hand the manifold H 2 x Si has two 
non compact directions; it is then meaningful to study the asymptotic behaviour of 
the Green's function at criticality, when the distance among points on H 2 diverges. 
We will also study such a behaviour in the limit of zero temperature (no compact 
directions at all). 

The 2-point Green's function of the operator A1/2 + m 2 is given by 

G(z,e,z',6')= jH dt h H? {t;z,z')h s i{t;B,e')exp{-m 2 t} . (8.1) 

Since we are interested in the asymptotic behaviour of the Green's function in the 
H 2 direction we can fix 6 = 6'. Moreover, at criticality it is m = 0. Using ( [7.91) 
and ( [7.1 2| ) we have then 



~, , ns rd 1 1 f°° , , /9 f 00 , w cosh(u>/2) 

G(z,z',6) = H / dtr 3/2 dw- 



(2tt) 5 / 2 v / cos h(d/ r ) + 1 ^0 Jd/r y/coshH - cosh(ci/r 
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~ f w 2 r 2 ( W t\ 
x £ex P -_-(„+- - . (8.2) 



n=0 



4t V 2/ (3 



After performing the integral in t and the sum over n we get 

G(s, 9) = P 1 r dw W eosh( W /2)cosech( W r/2/?) (g 3) 

4V27T 2 Jcosh(d/r) + 1 Wcosh(w) - cosh(d/r) 



In the limit d — ■> oo, with r and /3 finite, we can approximate the integral as 

-1/2 



u> 

r 



exp . (8.4) 



We finally get 



G mxS i(z,z',e)~ - -i expZ-jfi + lU . (8.5) 



As one can see, the correlation length at criticality is finite 

« = 2 (H)" ' (8 ' 6) 

namely, it is proportional to the two scales of the theory, (3 and r. 

To analyze the zero temperature case one cannot use the results (|8.5|) . 



which were obtained assuming (3 finite, but rather Eq. (|8.3| ) in the limit (3 — > oo. 
Thus we get 

g(».z'.«.)=4 / 1 r- , cmh(w/2) ■ M 

r7r v / cosh(c//r) + 1 Jd/r w yJcoBh.{w) - cosh(d/r) 

where x a = x' Q is the coordinate on R. We now take the limit d — > oo, r finite, and 
we get 

G H 2 xR (z,z',x ) -^=exp{-— 1 . (8.8) 

Tx 2 \]rd { 2r J 

As we can see, even in the case of a manifold which is non compact in all directions, 
we get for the asymptotic Green's function an exponential behaviour, and hence a 
finite correlation length, which is proportional to the radius of curvature r. 



9 Conclusions 

In this paper we have analyzed the large N limit of the 3-d Gross-Neveu model on 
the manifolds S 2 x Si and H 2 x Si. The physical observables can be regarded as 
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functions of two parameters, the coupling constant of the model and the inverse 
temperature (5. Thus the critical behaviour can be studied with respect to both. 
In particular, we have evaluated the free energy, the correlation length and the 
spontaneous magnetization at the UV critical value of the coupling constant, in 
the (^-function regularization scheme. These quantities completely characterize the 
thermodynamical properties of the model once the coupling constant has been 
fixed. For both manifolds, the numerical evaluation shows that the free energy 
density at the critical coupling is a smooth function of the temperature; hence no 
phase transition seems to occur as the temperature (1//3) varies. 

The asymptotic behaviour of the correlation function at the critical coupling is 
also obtained. We find, in complete analogy with the results of Ref.|2|, [| for the 
non linear a model, that the correlation length is finite even for manifolds which 
are non compact in all directions (H 2 x R). More precisely, the finite size effect 
is due to the non vanishing curvature of the manifold which introduces a length 
scale in the theory. Finally, we observe two major differences with respect to the 
non linear a model. First, the GN model appears not to be conformally invariant 
at criticality (in the large N limit). Then, we do not observe curvature induced 
symmetry breaking, that is, m c and b c stay zero even on curved background. 
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A The algebraic method 

- The space S 2 x Si - 



To solve Eq. (|6.4j) we apply the method described in Ref . ||26|| . According to this 
method we can construct the eigenvector of the Dirac operators by means of the 
eigenvectors, fa, for the scalar case 

&ofa = -V M 0"&, = ufa , (A.l) 

and the so-called covariantly constant spinors, ef , which are defined in this space 
by the following differential equations 



i 



V^±— 7/zJef = with p, = 1,2 (A.2) 
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V 3 e 



± 







(A.3) 



Note that, the above requirements on the spinors ef are compatible with the 
Bianchi identity which reads 



[V M) V,] ef = ^R^Yef . 



(A.4) 



The index i just labels the independent solutions of Eqs. (|A.2|) , (|A.3|) . As far as 
the equation ( |A.3j ) is concerned, it only says that the spinors ef do not depend on 
9. Then we only have to solve the equations (|A.2| ) on S%. The solutions of this 
problem have been obtained in Ref.[|27fl and are 



e i (t, X) = exp <^ — 



IT 



4( r >x) = ex p 



-IT 



sin(x/2 + tt/4) 

=F cos(x/2 + tt/4) ^ 
cos(x/2 + tt/4) 

. ±sin( X /2 + 7r/4) 



(A.5) 



(A.6) 



They satisfy the normalization condition 



(A.7) 



with i,j = 1,2. It is worth pointing out the right periodicity property of e 1( ^ 
spinors 

ef(T + 27i,x) = -et(T,x) • (A.8) 

However, the same does not occur for the variable 9. To solve this problem in the 
reconstruction method we use, for the scalars, twisted boundary conditions along 
9 



4> uj (t + 2ti,x,9) = (f> u (r,x,0) 
<!>u(t,x,0 + 2vr) = -<I> u (t,x,0) 



(A.9) 
(A.10) 



The eigenfunctions of ( |A.1|) with the above constraints are then 

<h™(T,x,9) = N lm exp{i(n+^e} X ,r) , (A.ll) 

where denote the spherical harmonics, n G Z, I = 0, 1, .., and —l<m<l and 
Ni m is a normalization constant. The corresponding eigenvalues are 

1 1 / 1\ 2 

ui n = ui + uj n = —l{l + l) + — \ [ n+-J . (A.12) 
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In terms of the expressions ( |A.5| ) , ( |A.6| ) and ( |A.11[ ) we can define four independent 
spinors on S 2 x Si, for fixed j, namely, 



|1) = 4>lmn4 > I 2 ) = ^(dfl&mn)^ , |3) = (plmn.tj , |4) = (dfi&mn) 



6: 



(A.13) 

where p, — 1, 2. It can be can be shown that the action of Dirac operator is closed 
on the space spanned by these four vectors, so that we can find combinations of 
them which are eigenvectors of the Dirac operator. In fact, the Dirac operator can 
be represented as a 4 x 4 matrix acting on the subspace spanned by the four vectors 
reported in Eq. (|A.13|) 



2n+l 



2/3 2 



V o 

The eigenvalues iX are then 

l^ 2 



1 


2n+l 




2/3 2 













2n+l 


1(1+1) 


2/3 





A 



n + 







1 | (l + lf 



\ 

2n+l 

2/3 

1 

o J 



1\ 2 1 I 2 

n+ 2) J 2+ 7 2 



(A. 14) 



(A.15) 



Note that there are no zero modes. Hence, the eigenvalues of the Dirac operator, 
zA^j, have the form 



\ 



1\ 2 1 (I + l) 2 
n + -) -7^ + ~ ^r- with n G Z, I = 0,1,.. 



2 



(A.16) 



It can be checked, writing the eigenvectors as linear combination of the basis vectors 
( |A.13| ), that the degeneracy is 2(1 + 1). 



- The space H 2 x Si - 



As for the space S 2 x Sp, also in the case of H 2 x Sp we look for the solutions of 
the equations 

1 



V, e 



V^±— 7p)ef = with fi = l,2 







(A. 17) 
(A.18) 



We find it is easier to solve the above problem in coordinates \ G R and r G [0, 27r] 
which are connected to x and y of ( |7.2j ) by 

sinhx sinr 1 



cosh x + sinh \ cos r 



cosh x + sinh \ cos r 



(A.19) 
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The independent solutions in this case result to be 

sinh(x/2) 



ef(r, X ) 



IT 

ex Piy 



=Fi cosh(x/2) 
cosh(x/2) 

e 2 ( r > X) = exp ( 

=pi sinh(x/2) 

The scalar problem yields to a continuous spectrum p3, E 



(A.20) 



(A.2i; 



with eigenf unctions 



and degeneracy 



U k + L0 Tl 



1 



1 



1 



n 



Pj 71 (cosh x) exp 



n 



-10 + imr 
2 1 



/x(fc) 



1 



7T 



9 IWfc - 



tanh ir\ Uf, 



(A.22) 



(A.23) 



(A.24) 



P; m (coshx) are the associated Legendre functions, with I = —1/2 + ik. 

As for S% x S^j, also in this case one can construct the four spinors of Eq. 
(|A.13[) and the action of the Dirac operator is closed on the corresponding subspace 
spanned by them. Hence, the eigenvectors of the Dirac operator can be expressed 
as a linear combination of them. We could in principle proceed in the same way as 
for x Si but we choose a simpler approach, using some known results. It can 
be easily seen that all the four independent spinors of Eq. ( A.lij| ) take the form 

|«) ~ \<t>n) ® \Xk) , (A.25) 



(n + |") 9 are eigenf unctions of the scalar Laplacian 



where \<p n ) = (2irf3)~ 1/2 exp 
on Si with twisted boundary conditions, while \xu) are spinors depending on H 1 
coordinates only. Since the eigenvectors of the Dirac operator are linear combina- 
tion of the four spinors ( |A.13| ), we can also say that each eigenvector is of the form 
( |A.25p . The action of the squared Dirac operator 

A 1/2 (tf r 2 x Sj) = A 1/2 (H 2 r ) + Ao(5j) (A.26) 

is then factorized in the sense that 



A 1/2 \ij n ,k) = A 1/2 |xfc) + A \(p n ) 



(A.27) 



and the heat kernel is just the product of the two heat kernels on H 2 and Si (cfr. 
Ref. H]). This is all what we need to write down the (^-function. 
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B Poisson sum formula for S* 2 



The general expression for the Poisson sum formula (see Ref.EH]) is 



exp { + r]) 2 + 2m- (I + rj) 



l=-oc k 7 7 

r (x + / 7 ) 2 



7 



'Alrt 



J2 ex p 



l=—oo 



At 



— Iixilr] 



We put rj = and 7 = 27r and derive both the sides of (|B.1|) with respect to x; we 
get 



% 

2^ 



/ exp j— Z 2 t + z/xj 



1 f> +2 W> +M)2 



Z=— 00 

Then we observe that 



2tV47rt 



Z=— 00 



so that 



00 00 / 

^2 sign(m)e* mx = 1 + 2i sin(mx) = 1 + i cot f 

^=— 00 m=l 

sign(m) = — [ dx cot f-) e " ?mx . 
27T Jo \2/ 



4t 



(B.2) 



(B.3) 



(B.4) 



Now we multiply (|B.2|) by ^- cot(x/2) and integrate 
1 z^ 71 



4vr 2 



/" /x\ r i 

y da; I cot ^— J exp |— l 2 t + i/xj 



2 

(47Tt) 



2vr ^ 

j J dx cot f — J x (x + 27rZ) exp 
" i=— 00 



^x + 2?r/) ! 



4t 



(B.5) 



Then we use ( |B.4| ) and we get 

— ^ / sign(-/)exp{-/ 2 t} 



Z=— 00 



—r fdx cot f J) £ (x + 2vrZ) exp J - 



(x + 2%lf 



(B.6) 

(47rt)S J V2/ z — ' * [ 4t 1 

Performing the change of variable x' = x + 27r/ and exchanging the sum with the 
integral we have 

(x + 2tt/) 2 1 



f 2n f x \ °° f 

/ dx cot I — ) Yj (x + 27r/) exp < - 

J ° ^ 2J l=-oo { 



At 



dx x cot ( — 1 exp < > . 

\2J P l At I 



(B.7) 



22 



Substituting this expression in ( |B.6| ) we finally get the Poisson sum formula for the 
spin 1/2 heat kernel on S 2 

-£z exp{-/ 2 t} = -i-3- [°° dxx cot(|)exp(-^-} . (B.8) 
7r ^ (47rf)2 J-oo \2J { At J 
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Figure 1: S% x Sp. The free energy density, w c , of (|6.19|) is plotted as a function of 
P/r, up to the factor l/(167r 2 r 3 ). 




Figure 2: x Sp. The free energy density, w c , of (|7.22|) is plotted as a function 
of {3/r, up to the factor l/(167r 2 r 3 ). 
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